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We have performed ab-initio calculations of exchange couplings in the layered vanadates CaV^Os, 
MgV 2 05, CaV307 and CaV-iOg. The uniform susceptibility of the Heisenberg model with these 
exchange couplings is then calculated by quantum Monte Carlo method; it agrees well with the 
experimental measurements. Based on our results we naturally explain the unusual magnetic prop- 
erties of these materials, especially the huge difference in spin gap between CaV^Os and MgVjOs, 
the unusual long range order in CaVaOy and the "plaquette resonating valence bond (RVB)" spin 
gap in CaV^Ot). 
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The series of insulating vanadates with the compo- 
sition MV„02n+i is of considerable interest because it 
shows a variety of strange magnetic phenomena which 
are not completely understood. Among these quasi-two 
dimensional layered materials has originally at- 

tracted a lot of attention as the first two-dimensional 
material with a gap in the spin excitation spectrum 
attributed to a plaquette RVB state j^j . The other mem- 
bers of this family are also of interest. CaV~307 shows 
long range Neel order with an unusual ordering g] that 
could not be explained by simple models CaV^Os has 
a large spin gap of about 600 K ||, but for MgV 2 Os on 
the other hand, which has nearly the same lattice struc- 
ture, there is only contradicting evidence for a possibly 
gapped state with a much smaller gap of only about 20 
K §. 

Despite extensive theoretical work on these materials 
@,|J|Q their magnetic properties are not yet fully un- 
derstood. In contrast to planar cuprates, where a hole 
in x 2 — y 2 -orbitals of Cu results in a strong antiferro- 
magnetic exchange coupling for 180° bonds and a weak 
ferromagnetic one for 90° bonds, the interactions in these 
vanadates are much more complicated. Not even the sign 
of many of the exchange couplings is obvious in these ma- 
terials. Only straightforward numerical ab-initio calcu- 
lations can give information about relative and absolute 
values of the exchange couplings in these systems. 

In this Letter we report on such an ab-initio calculation 
using the LDA+U method [jl3| to compute the electronic 
structure and from it the exchange couplings. We then 
use the quantum Monte Carlo (QMC) method to calcu- 
late the uniform susceptibility assuming these exchange 
couplings, and compare it with experimental measure- 
ments. We find good agreement and are thus confident 
that our explanation of the magnetic properties of these 
materials captures the relevant physics. 

The LDA+U method was shown to give good results 
for insulating transition metal oxides with a partially 
filled d-shell The exchange interaction parameters 



can be calculated using a procedure based on the Greens 
function method which was developed by A.I. Lichten- 
stein pq , ^6| . This method was successfully applied to 
calculate the exchange couplings in KCUF3 |16| and in 
layered cuprates [jl7| . 

The LDA+U method is essentially the Local 

Density Approximation (LDA) modified by a potential 
correction restoring a proper description of the Coulomb 
interaction between localized d-electrons of transition 
metal ions. This is written in the form of a projection 
operator: 
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where | inlma) (i denotes the site, n the main quantum 
number, I- orbital quantum number, m- magnetic num- 
ber and a- spin index) are d-orbitals of transition metal 
ions. The density matrix is defined by: 



ImGf nlmMlm ,(E)dE, (2) 



where G a . , ,(E) 



(inlma | (E — H) 1 | inlm a) 
are the elements of the Green function matrix, N a = 
Tr(n a mml ) and N = N 1 * + N l . U and J are screened 
Coulomb and exchange parameters. The U mm 'm"m'" 
is the screened Coulomb interaction among the nl elec- 
trons which can be expressed via integrals over complex 
spherical harmonics and U and J parameters. For the 
CaV„02n+i compounds the values of these parameters 
were calculated to be £7=3.6 eV, J=0.88 eV via the so- 
called " supercell" procedure |t8| . The calculation scheme 
was realized in the framework of the Linear Muffin-Tin 
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Orbitals (LMTO) method jL9| based on the Stuttgart 
TBLMTO-47 computer code. 

The inter-site exchange couplings were calculated with 
a formula which was derived using the Green function 
method as second derivative of the ground state en- 
ergy with respect to the magnetic moment rotation angle 
II® 

T - = P ,V j P 

"V / j mm' Amm'm"m"' m"m"' \"J 
{m} 

where the spin-dependent potentials / are expressed in 
terms of the potentials of Eq. ([!]), 

P , = V lX , - V il , (4) 

-*mm' v mm' 'mm' \ J 

The effective inter-sublattice susceptibilities are defined 
in terms of the LDA+U eigenfunctions ip as 
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The main building block of the crystal structures of the 
CaV„02n+i series is a V ion roughly in the center of a 
pyramid of oxygen ions. The lowest energy orbital is the 
V3d-orbital of xy-symmetry (using a convention where 
the axes of the coordinates system are directed toward 
the oxygen ions in the plane), which is the orbital whose 
lobes look at the directions, where the overlap with the 
oxygen is the smallest. The crystal field splitting between 
xy-orbital and other 3d-orbitals is so strong, that in the 
LDA band structure the xy-band is separated from the 
rest of V3<i-bands by a small energy gap jl(J . 




FIG. 1. The occupied d-orbitals of V 4+ ions in CaVsC^. 
Oxygen atoms (denoted as small black circles) form pyramids 
with V atoms inside them. 



As a consequence the degeneracy of the V3d-shell is 
lifted and the single d electron of V 4+ ion occupies this 
xy-orbital, which reminds us of the cuprates, with a sin- 
gle hole in the x 2 — y 2 -orbital. The important difference 
is that while in cuprates all copper atoms are in the same 
(x, y)-plane as the x 2 — y 2 -orbital, in these vanadates the 
vertices of the pyramids point alternatingly up and down 
with respect to the basal plane. Thus the V ions in their 
centers are correspondingly above and below the central 



plane, as can be seen in Fig. [i]. As the xy-orbitals are par- 
allel to this plane, the overlap (and hence the exchange 
couplings) are expected be stronger for vanadium ions 
situated on the same side of the plane. We will show 
that this is indeed the case. In addition to this alterna- 
tion, a tilting of the pyramids is present in the crystal 
structure of these compounds, which seriously influences 
the interactions. 

Another important difference with the cuprates is that 
the a:?/-orbital has a 7r-overlap with the in-plane oxygen 
atoms in contrast to a much stronger c-overlap in case 
of Cu 2+ . Consequently one can expect much weaker ex- 
change interaction in vanadates as compared to cuprates. 
All the more surprising is thus the fact that the spin gap 
in CaV205 (616K Q), is larger than the typical values 
for the similar cuprate ladders (« 460 K [po| ). 

Equation (|J) was derived as a second derivative of the 
total energy with respect to the angle between spin di- 
rections of the LDA+U solution. The LDA+U method 
is the analogue of the Hartree-Fock (mean-field) approx- 
imation for a degenerate Hubbard model |0. While in 
the multi-orbital case a mean-field approximation gives 
reasonably good estimates for the total energy, for the 
non-degenerate Hubbard model it is known to underes- 
timate the triplet-singlet energy difference (and thus the 
value of the effective exchange coupling ) by a factor of 
two for a two-site problem (Ehf = ^j" and E exact = 
where t <C U is inter-site hopping parameter). In the 
problem under consideration there are two types of con- 
tributions to the exchange interaction parameters Jy. 
The first one is due to the xy — xy orbitals hopping, and 
as only this orbital is half-filled this contribution directly 
corresponds to the non-degenerate Hubbard model and 
its value must be multiplied by a factor of two to correct 
the Hartree-Fock value. Other contributions are due to 
the hoppings to all other orbitals and as the mean-field 
approximation is much better for multi-orbital model this 
part can be used unmodified. 



TABLE I. Calculated exchange coupling parameters (K). 
"Minus" sign means ferromagnetic exchange. 
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In the crystal structure of these MV„02n+i com- 
pounds layers are formed by the VO5 pyramids sharing 
edges of their bases. The V atoms of the layer form 
roughly square lattice which is 1/ (n+l)-depleted (Fig. ||). 
As mentioned earlier, the strongest interaction must be 
between V atoms which are situated on the same side of 
the plane (above or below) and are marked by the light 
gray and the dark gray circles in Fig. |[ 
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FIG. 2. The basic crystal structure and the notation of exchange couplings in CaV20s and MgV20s (left panel), CaVs07 
(middle panel) and CaV4C>9 (right panel). V atoms represented by large circles with different colors have different z-coordinate. 
Oxygen atoms are shown by small circles. The long range magnetic structure of CaW^Or is depicted by white arrows. 



In the MV2O5 crystal structure these atoms form lad- 
ders with interactions along the rung and the leg of the 
ladder denoted as J2 and J3 and interaction between lad- 
ders as Ji (the notations are chosen to reflect the inter- 
atomic distances; the shortest one is between atoms on 
different sides of the plane). In CaVaOy the atoms on the 
same side of the plane form zigzag chains with interaction 
inside this chain denoted as J3 and interactions between 
chains as J\ and J 2- For CaA^Og such atoms form "meta- 
plaquettes" (this term was proposed by W.Pickett jn| to 
distinguish them from the plaquettes formed by atoms 
with the shortest V-V distance). The interaction inside 
the metaplaquette is denoted by J3, between metapla- 
quettes on the same side of the plane as J4, and between 
metaplaquettes on the opposite sides of the plane by Ji 
and J2- 

Our calculated values of the exchange couplings are 
presented in the Table I. It can immediately be seen 
that indeed the strongest interactions are between atoms 
on the same side of the plane ( the ladder exchanges 
Ji, J3 for CaV205 and MgV20s, the zigzag exchange 
J 3 for CaV3C>7 and the metaplaquette exchange J3 for 
CaV40g). However all other parameters are not negligi- 
ble and for CaA^Og the frustrating interactions between 
metaplaquettes are not less than 60% of the value of in- 
teraction inside the metaplaquette. 

For a comparison of these exchange couplings to mea- 
surements on the materials and to discuss the magnetic 
properties we consider the temperature dependence of 
the uniform susceptibility x{T)- This quantity depends 
very sensitively on the exchange constants and can be 
both easily measured in experiments and calculated nu- 
merically for the Heisenberg model. To calculate the 
x(T) we use the continuous time version of the quantum 
Monte Carlo (QMC) loop algorithm ^ . This algorithm 
uses no discretization of the imaginary time direction and 
the only source of systematic errors are thus finite size 
effects. The lattice sizes were chosen large enough so that 
these errors are much smaller than the statistical errors of 
the QMC simulations. All of these simulations, with the 
exception of the model for CaVsOy suffer from the neg- 



ative sign problem due to frustration effects. Improved 
estimators |2^] can be used to lessen this sign problem, as 
was done in Ref. |l2| for similar models. The frustration 
is especially strong in MgV 2 05 and CaV/tOg. In these 
two cases x(T) could not be calculated to temperatures 
far below the susceptibility maximum. However x(T) 
at intermediate temperatures is already very sensitive to 
the values of the exchange couplings and a comparison is 
thus possible. 
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FIG. 3. Comparison of the temperature dependence of the 
uniform susceptibility \{T) measured in experiments to that 
calculated in QMC using the LDA+U values for the exchange 
coupling of a) CaV 2 5 , b) MgV 2 5 , c) CaV 3 7 and d) 
CaV409. The ^-factor of V +4 was taken to be g = 1.96 



The g-factors for V 4+ (d 1 ,S'=l/2) compounds are ob- 
served to be in a narrow range about g — 1.96 [ p4| . We 
use this value to compare (in Fig. ||) the calculated x(T) 
to measurements done by Isobe and Ueda |23| over a 
wide temperature range. Taking note that the purity of 
the samples, and thus the normalization of the experi- 
mental data, is not known precisely and keeping in mind 
how sensitive is x(T) dependence to the values of the 
exchange couplings we find reasonable agreement of the 
calculated x(^) with the experiments. 

The agreement is especially good in the case of 



3 



Cav^Os where the exchange coupling J 2 on the rung 
of the ladder is dominant. Our LDA+U results thus con- 
firm the weakly coupled dimer picture proposed for this 
compound in 

The magnitude of the exchange coupling however de- 
pends strongly on the tilting of the oxygen pyramid. The 
compound MgV^Os has nearly the same crystal struc- 
ture as CaV2C>5, but as the Mg ion has a smaller ionic 
radius than Ca ion, the tilting of the oxygen pyramids in 
MgV^Os is stronger. This causes the exchange coupling 
on the rung Ji to be more than four times smaller and ex- 
plains the huge difference in the uniform susceptibilities 
and spin gaps of these two compounds. 

The exchange couplings in MgV205 are all of the same 
order, which positions this material in the strongly frus- 
trated region of the trellis lattice phase diagram p5| , pre- 
cisely the region about which not much is known yet. 
Even the experimental evidence is ambiguous, with dif- 
ferent claims about the existence or non-existence of a 
small spin gap reported in the literature More de- 
tailed experimental investigations are in progress, which 
will help to understand not only this compound but also 
the properties of the strongly frustrated trellis lattice 
Heisenberg model. 

Of the compounds under consideration only CaVa07 
has no frustration, and at low temperatures our QMC 
simulations gave a magnetically ordered Neel state, with 
the magnetic structure coinciding with the one observed 
in experiments. This structure, with ferromagnetic or- 
der on the short 3-atom "rungs" and antifcrromagnetic 
order along the infinite chains could not be easily ex- 
plained previously. The most surprising result is that 
the exchange couplings Ji and J± which were thought to 
be equivalent Q] are actually of opposite sign, most prob- 
ably due to the tilting in the lattice structures (Fig. 1). 
These exchange couplings also naturally lead to the ob- 
served experimental order, and no recourse to quantum 
fluctuation effects, as proposed by Kontani et al. Q, is 
necessary. While there is qualitative agreement between 
the calculated and measured susceptibilities it is not as 
good as in the other compounds, due to closeness to a 
ferromagnetic state and large cancellation effects in the 
effective exchange couplings that increase the errors of 
the LDA+U exchange coupling estimates. 

Of the four compounds CaV^Og is the most stud- 
ied one. We find that indeed the originally proposed 
plaquette-RVB state p) is relevant, but with the larger 
meta-plaquettes (J3) being dominant, as suggested by 
Pickett |l(| and by Kodama et al || . However, in con- 
trast to the estimates of [|| we find that the coupling 
J4, while smaller than J3, is comparable to the other 
exchange couplings. 

Our estimated exchange couplings can explain the or- 
der of magnitude of the spin gap of CaV^g. With our 
values of J's the unfrustrated J3-J4 model was shown 
to have a spin gap 0] of about 60K, which is further 



increased by frustration due to J\ and J 2 . While ex- 
act theoretical estimates for the spin gap in this strongly 
frustrated regime are not available, this estimate agrees 
with the measured spin gap of f 07K, and the good agree- 
ment between the measured and predicted uniform sus- 
ceptibilities is a further indication for the validity of our 
estimates. 

Additional comparisons can be done by calculating the 
magnon dispersion of the model and comparing it to ex- 
periments, as was recently done by Mambrini and Mila 
PJ. According to their data, obtained on a small f6-site 
cluster our ratio J3/J4 might just be slightly too large 
to give the dispersion minimum at (0, 0), but further cal- 
culations on larger clusters are necessary to clarify this 
point. 

To summarize, using LDA+U calculations of the ex- 
change couplings for a series of layered vanadate com- 
pounds we can explain the puzzling magnetic properties 
of these materials and find good agreement between the 
ab-initio predictions for the uniform susceptibility and 
experimental measurements. Such calculations are espe- 
cially needed for systems such as the vanadates inves- 
tigated here, where not only relative strength, but even 
the sign of the exchange interaction can not be estimated 
using simple geometrical arguments. Applied to a series 
of layered vanadates we can explain the difference in spin 
gap of CaV205 and MgV^Os, obtain a natural explana- 
tion of the unusual magnetic order in CaVs07 and con- 
firm the plaquette-RVB state of CaV40g with ab-initio 
calculations. 
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